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In this paper we study the question of effective field assignment to measured or nonperturbatively 
calculated spectral functions. The straightforward procedure is to approximate it by a sum of 
independent Breit-Wigner resonances, and assign an independent field to each of these resonances. 
The problem with this idea is that it introduces new conserved quantities in the free model (the new 
particle numbers), therefore it changes the symmetry of the system. We avoid this inconsistency 
by representing each quantum channel with a single effective field, no matter how complicated 
the spectral function is. Thermodynamical characterization of the system will be computed with 
this representation method, and its relation to the independent resonance approximation will be 
' discussed. 
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A fundamental field theory is defined through its microscopic field content, representing the fundamental degrees of 
freedom (dof), and the interaction between them. But microscopic dof are not always observable: in low energy QCD 
we find hadrons instead of quarks and gluons. Also in an ordinary gas we might encounter rather molecules than 
I separate electrons and nuclei (even nucleons). Although this statement reflects a trivial fact, it is hard to consistently 
Oh> implement the rather smooth shift/overlap in the nature of the relevant dof in an interacting theory. If we want to 
' work directly with the observable dof, we need to build effective theories. 

A primary principle of construction of the effective theories is that the symmetries of the fundamental system should 
be present unchanged in the effective model. This ensures that they belong to the same universality class, ie. they 
describe the same physics. In fact, the by-now textbook examples of the QCD effective models like the nonlinear or 

■ linear sigma models, chiral sigma models, Polyakov loop models, Walecka model, etc, use this strict guiding principle 
| to represent low-energy forces generated by QCD 0, Q ■ In these models one defines the field content to represent the 

■ low-lying bound states of QCD, and construct a Lagrangian which respects the symmetries of QCD. The parameters 
' of the Lagrangian are eventually tuned to reproduce most observed features of the particle spectra. 

In the examples listed above effective fields represent not just the (lowest lying) observable hadrons with the 
corresponding quantum numbers but actually the complete quantum channel. We get into trouble, however, when 
^sO . we want to represent individually also the higher energy states in a given channel; in case of QCD the plethora of 
excited hadrons. The representation of the observed particles would require to assign a field to each hadron. The 
manifestation of this idea, treated in the ideal gas approximation is the widely used hadron resonance gas (HRG) 
model, which indeed provides surprisingly accurate description to certain quantities in QCD in the low temperature 
regime For a more realistic description^ in particular to describe transport properties, one has to consider also 

the widths of the peaks (off-shell transport) 0, H[ • The problem of conservation of energy-momentum tensor and other 
conserved currents can be consistently solved in the framework of Kadanoff-Baym equations [9]. In this extended 
version one still assigns a field (Breit-Wigner resonance) to each peak of the measured spectral function separately 

In these approaches, however, one gives up the principle of keeping the symmetries of the fundamental theory, 
by generating a large number of new conserved quantities at tree level, namely the particle numbers of the excited 
hadrons. With other words: each quantum channel is represented multiple times. If we do not want this unphysical 
extra symmetry to appear on the level of observables, the tree level contributions must be compensated by the effect of 
interactions. This may require, however, strong radiative corrections of 0(1), eventually resulting in a nonperturbative 
system. Then tree level results cannot be considered reliable, whatsoever nice is the arising description. This situation 
would question the effective nature of the chosen model. 

There are cases when the idea that we assign a field (a new dof) to each peak of the spectral functions indeed 
leads to false conclusions. A well-known example is the hydro gen atom: if we took all bound states into account as 
independent dof then we would obtain false thermodynamics [ill . QjJ ■ In the standard solution [ll], QjJ we quit at 
some scale the framework of quantum mechanics, and explicitly assume that bound states with certain size cannot 
exist in the matter. Although physically motivated, nonetheless this requirement is not self-consistent, and it is hard 
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to imagine that it can be phrased formally on the level of the Lagrangian. In a similar spirit, one can try to improve 
HRG by cutting off the number of hadronic resonances above an appropriately chosen energy scale (l3j . 

We can study this question within the scattering theory, too. In the quantum mechanical approach, bound states 
appear as a jump in the scattering phase. Their contribution to free energy is described by Beth-Uhlenbeck formula 
}12| provided the scattering phases of resonances are simply additive |14j . In the generic situation, however, this is not 
the case. Finite width peaks should be included in the S-matrix as Breit-Wigner resonances with complex amplitudes 
[TH fl6|. The unitarity of the S-matrix then yields a consistency constraint, which should be solved. The resulting 
amplitudes are thus results of interference effects which might decrease the "weight" of the individual finite lifetime 



dof. If the finite width peaks of the spectral function are well separated, then a careful analysis can show 17| that 
the additivity of their scattering phases is indeed true. 

In this paper we approach this question from the field theory point of view. The potential problem coming from 
the introduction of new conserved quantities can be avoided by simply representing each quantum channel by a single 
field - no matter, how many peaks are in the corresponding spectral function. As we have discussed above, associating 
a single field with the singled out quantum channel works well in a lot of effective models of QCD. The representation 
of a complicated spectral function by a single field, technically is similar to the perturbative 2PI technique [lH , where 
the propagator (e.g. the spectral function) of a channel is a single dynamical object. However, in this case we go 
beyond the reorganization of the perturbative series, we shall construct a non-perturbative scheme where the spectral 
function is the input. We expect (and actually will prove) that our scheme shares with the 2PI technique a number 
of good features like energy-momentum conservation, causality, unitarity and renormalizability. 

In this paper we try to use this idea to build up a consistent one-dof bosonic effective model. We start with a 
very simple analysis, demonstrating that the correct representation of the spectral function is relevant (Section [TTJ) . In 
Section Mil we present the one-dof model and prove its consistency. In Section IIVI we calculate the thermodynamical 
observables (energy density, pressure etc.) from this model, and compare the results with the ones assuming inde- 
pendent excitations. The comparison is facilitated by a prescription how to measure the "number of dof" . We will 
see that overlapping dof are indeed not independent, and they always behave as if they were less than one degree of 
freedom. In Section PVTl we give the conclusions. 



II. WHY THE REPRESENTATION OF THE SPECTRAL FUNCTION IS IMPORTANT? 

The goal of this Section is to demonstrate that the spectral function alone does not determine the physics completely: 
the symmetry properties influence certain observables, in particular the thermodynamics crucially. 

First let us define the spectral function: we choose an operator A with fixed quantum numbers, and compute 
the QAAf(t) — (l/K*)) -^(Q)]) expectation value, and take its Fourier transform Qaa^(po)- If the system is spatially 
translation invariant then the spatial momentum is part of the quantum numbers, too. As it is well known [19f . this 
expectation value at zero temperature is proportional to the density of states, at finite dimensional quantum mechanical 
systems it consists of discrete Dirac-deltas at the position of the energy levels. Choosing different operators to measure 
the spectrum results in different weights of the Dirac-deltas, but their position will be the same. 

Now consider a one-dof free scalar Klein-Gordon model. Its Lagrangian reads 

C = \<PK<P, K = -d 2 t -u 2 k , ^ = k 2 + m 2 . (1) 

We can use the operator tp(t, k) to measure the spectral function in the momentum k, spin quantum channel. It 
reads j?(po) = 27r sgn(j0p)(5(p n — ufj). We remark that because of the canonical commutation relations this spectral 
function satisfies a sum rule [19j . 

As a second example we take two scalar degrees of freedom ip\ and ip 2 with the same mass, and construct their free 
model: 

r-W*. (2 , 

To measure the spectral function in the momentum k, spin quantum channel, we can choose for example the operator 
A{t, k) = ipi(t, k) cos 9 + (f2(t, k) sin# for arbitrary angles. The measured spectral function will be, independently on 
the angle, Qaa*(Po) — 27r sgn(p )<5(p 2 — w fc)> the same result as before. 

In this model there is an extra internal symmetry, the 0(2) = U(l) rotation of the two fields. If this symmetry is 
physical, then there may exist a measurable quantity which makes a distinction between the two fields; sometimes it 
is the electric charge attached to the U(l) rotation. But if we introduce the two fields exclusively to represent the 
spectrum, the extra symmetry is non-physical, and the physical operators do not see the difference between the two 
degrees of freedom. In this case physically there exists only one quantum channel with k momentum and spin 0. 
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Now we have two models which give the same spectral function. The thermodynamics, described by these models, 
is still different. For example the Stefan-Boltzmann (SB) limit for the energy density would be e = 7r 2 T 4 /30 in the 
first case and twice as much in the second case. This can be physical, of course, if the U(l) symmetry has a physical 
background. But if the two degrees of freedom were introduced exclusively to represent the spectrum, then the extra 
symmetry is nonphysical, and we have false results for thermodynamics, despite the fact that we correctly described 
the spectrum. 

For practical reasons we may still want to choose a two-particle representation even if the physics do not support 
the new symmetry. For example if we have two Dirac-delta peaks in the spectral function then, according to (171 ] , we 
may use a two-particle representation. In the limiting case when the two peaks merge, we arrive to the second model. 
Since the nonphysical symmetries appear in the observables, in particular in thcrmodynamical quantities, to a correct 
two-field description we have to choose an interaction which breaks the artificial symmetry. Now we consider two 
examples for this symmetry breaking interaction within the framework of a quadratic model, where the interaction 
appears as off-diagonal elements of the kernel matrix. 

Consider for example the following model: 



This Lagrangian has no more U(l) symmetry. A special property of this model is that the kernel has no inverse. 
This means that one degree of freedom (now the ipi — if 2 mode) has no dynamics, it completely decouples from the 
system. The remaining mode if\ + if 2 has the same dynamics as in (JTJ), so the physics is the same too. In particular 
the spectral function (up to a normalization facto^) and the thermodynamics is the same. 
Another example is the model 



The two eigenvalues are K and K — 2A, so the spectral function exhibits two peaks. But in case A — > 00, the mass of 
one eigenstate remains m, while the other has a mass-squared m 2 + 2A — > 00 which decouples from the system. So 
the measurable part of the spectral function as well as the thermodynamics at any finite temperature behave as if we 
had one dof. 

So apparently we can work with more dof than the physical ones. The price we have to pay for it is that we have 
strong interactions. This must be so, since in the non-interacting (diagonal) case the SB limit is twice as much as 
in the interacting case: so the interactions give 0(1) corrections relative to the free case. The lesson is that wrong 
assignment of the dof 's alone may lead to an apparently strongly interacting theory, even if with correctly chosen 
fields the theory is free. 



III. REPRESENTATION OF A GENERIC SPECTRAL FUNCTION WITH A SINGLE FIELD 

Let us try to generalize the observations of the previous Section for the case when we observe, after fixing all possible 
quantum numbers, a generic spectral function. Let us assume that we can identify finite width peaks. Heuristically, 
the width of the peaks can be treated as a "resolution power". If, according to this resolution, the separation of 
the peaks is large, then earlier studies [ill E3 tell us that these behave as distinct objects: the independent particle 
approximation works well. If the separation of the peaks starts to be comparable with their width, then by their own 
resolution they are more and more like one object, and we arrive at the case discussed above. Then the lesson of the 
examples of Section [IT] tell us that if we further insist to a multi-dof description, the theory which correctly describes 
physics becomes strongly interacting, nonperturbative. 

This complication can be avoided if we do not introduce new dof into the system, and assign a single field to each 
quantum channel: then the symmetries of the original system is not altered. As a price, however, we have to work with 
a more complicated quadratic theory which can represent a generic spectral function. In this Section we construct 
a 1-dof system which can represent a generic spectral function, then we demonstrate that it describes a consistent 
theory, if the input spectral function comes from a physical system. 



1 We can determine the spectral function measured by the the above defined A, it results in qaa{po) = 7rsin 2 (0 + J) sgn(po)<5(pg — ^l)- 
This depends on 8 since the action is not U(l) symmetric 




(3) 




(4) 
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A. The Lagrangian 



Let us assume that we have measured a scalar spectral function, and there are no other hidden quantum numbers 
in the system. This means that the system must be represented with a single scalar field ip. We are looking for a 
quadratic Lagrangian which could represent the given spectral function (interactions can be necessary to correctly 
describe spectral functions of other quantum channels, which is not the topic of the present paper). We want moreover 
a system where energy and momentum are conserved, so we want space and time translation symmetry. Then the 
most generic scalar field theory can be written as: 

S[ip] = ^J d A xip{x)JC{id)ip{x). (5) 

In Fourier space we can write 

sM = lj-0f<P*(pMpMp)- (6) 

If K, is not a polynomial, then this action is not local. Since ip(x) is real, (p(—p) = (p*(p); since <S* is also real, the 
kernel must satisfy JC*(p) = lC{—p) = K,{p). 
The spectral function is defined as 

Q{x- y) = ([ip(x),<f(y)]) . (7) 

Since the theory is quadratic, so the commutator [ip(x) , (p(y)] is a c-number, thus we can also omit the expectation 
value. Using the standard procedure [19j . we compute the retarded retarded Green's function (Landau prescription) 
as 



G R (p) = Kr 1 (j> + ie,p). (8) 
The spectral function is the discontinuity of the retarded Green's function 

g(jp) = DiscGij(p) = lim [/C _1 (po + ie,p) - JC^ 1 {p - ie,p)] . (9) 

That means that a kernel determines the spectral function in a unique way. 

In our approach the spectral function is the input, and we want to construct the corresponding kernel. The retarded 
Green's function can be obtained via the Kramers-Kronig relation: 

G R (p)=[p^±-. (10) 
J 2tt pa - lu + is 

Using ©, the kernel can be determined as 

'p^lY 1 , (id 

where P means principal value. That means that the kernel is also completely determined by the spectral function. 



B. Consistency 

When one works with a nonlocal theory like ([SJ, one has to worry about the consistency of the model which, in 
case of a field theory, means that the requirements of energy- and momentum conservation, causality, unitarity, and 
Lorentz-invariance (in case of zero temperature) must be satisfied. 

First we remark that within perturbation theory it is also possible to change the tree level propagator and work 
with a self-consistently determined generic propagator instead: this is the 2PI resummation [18|]. This method is 
known to be a fully consistent way of treating the system perturbatively. Since one can formulate 2PI resummation 
as a model with nonlocal kernel [20j, the present method can be thought to be a nonperturbative realization of the 
consistent 2PI resummation method. 

But, independently on the perturbative treatment we can prove the consistency criteria one-by-one. 
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a. Energy and momentum conservation This requirement is satisfied if the action is time and space translation 
invariant. From the definition ([5]) this fact is trivial. 

b. Causality: Causality means in relativistic invariant systems that measurements at space-like separated points 
must not influence each other (in non-relativistic systems causality is not an issue) . Mathematically this requires that 
local operators at space-like separated positions must commute. We show that if the input spectral function is causal 
then the system is causal, too. 

A local operator built up from a single scalar field has the form 0(x) = Yii—i Ai, where Ai = a,i(x,id)ip(x) where 
a,i(x, id) is a local differential operator (ie. at most polynomial in the derivatives). We define Q(y) in a similar way 
as YiiLi Bi- Then we can use the following identity 

n m 

[6{x),Q{y)\ = ££.Bi • • • B j _ 1 A 1 . . . Ai.tlAi, B 3 ]A l+1 . . . A n B J+1 ...B m . (12) 
i=i j=i 

The commutator can be expressed as [j4j,Bj] = cij(x, id x )bj(y, id y )[<p(x), <p{y)]- Since aj and bj are local differential 
operators, an infinitesimally small opened neighborhood of the point x is enough to compute them from the field (p. 
But if x and y are spatially separated points, ie. (x — y) 2 < 0, then a sufficiently small opened neighborhood is also 
spatially separated. Since g(x — y) = [tp(x),<p(y)] is the input spectral function, it is guaranteed that it is zero for 
spatial separation. This means that in a small neighborhood it is identically zero, therefore all derivatives are zero. 
This means that is zero Vz,j, and finally [0(x),Q(y)] = 0, too. 

c. Unitarity Physically unitarity means that there is no loss of probability in the system. If we use, for example, 
finite lifetime particles, then the time evolution will not be unitary, since the decaying state will be missing from the 
final states. To circumvent this problem we have to keep all possible energy states, then a decaying state turns into 
another state and the complete probability will be the same. As a dynamical quantity, therefore, we have to work 
with the complete spectral function. 

Mathematically one can examine the unitarity (and the corresponding reflection positivity for euclidean systems) 
of models where the kernel is an arbitrary polynomial [2l| . The result is, formulating in the language of the spectral 
function, that for unitarity the spectral function must be positive for all positive frequencies q{uj > 0,p) > 0. Since 
we use the spectral function as an input stemming from a real quantum field theory, this requirement will be satisfied. 
Therefore unitarity is also granted. 

d. Lorentz invariance From the definition ([5]) it is clear that if the kernel is Lorentz invariant, then the action 
is Lorentz invariant, too. For the reverse argumentation, in Lorentz invariant systems the input spectral function 
leads to a Lorentz invariant retarded Green's function, which leads to a Lorentz invariant kernel. So, for a physically 
sensible zero temperature spectral function the Lorentz invariance of our model is also granted. 

The input spectral function, if it comes from a finite temperature system, need not be Lorentz invariant. Still, in 
most part of this paper we will use a Lorentz invariant input. 

We see, that all requirement of a consistent field theory could be satisfied with the Ansatz ([5]), provided we use 
an input spectral function which is physically sensible. Now we are in the position that we can use it for calculating 
physical observables. 



IV. THERMODYNAMICS 

Given the Lagrangian of the system, in principle we can calculate the expectation value of any physical observable, 
and, in particular, approach thermodynamics. But, because of time-nonlocality, some standard techniques do not 
work in this case. First of all we do not have canonical formalism (more precisely we would need infinitely many 
fields to represent the canonical formalism). Therefore there is not a direct connection between the imaginary time 
formalism and the thermodynamics: in fact, the standard derivation is based on the fact that the Hamiltonian density 
depends on the canonical momentum as n 2 /2 [l9| . Therefore thermodynamics should be started from something 
which is meaningful also microscopically: this can be the energy density. Known for all temperatures, this is then 
enough to calculate all the other thermodynamical quantities like pressure or entropy density. 



A. Energy-momentum tensor 

The energy density is the 00 component of the energy-momentum tensor, which represents the Noether currents 
belonging to the time and space translation invariance, respectively. We can use the standard procedure to determine 
the energy-momentum tensor: we consider position dependent translation a^(x) with which = x M + a^(x) and the 
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shifted field satisfies <p'(x') — <p(x). Then we can calculate the action with the shifted field S[<p'], for which we expect 
the behavior S[ip'] — S[<p] = J d xT^ v {x) dndv, where T^ v is the energy momentum tensor. 
Statement: The energy momentum tensor belonging to the action §5§ is 



Tp V {x) = -<f{x) D^K.(id)ip(x) 



where 



Da V IC(id) 



~dlC{p) 




dpi 1 


p—tid 



and the symmetrized derivative is defined as 

n 

f(x)[(idr] sym9 (x) = —j2[(-i d yf(x)}[(idr- a g(x)}. 



(13) 



(14) 



(15) 



a=0 



(16) 



Proof: We represent the kernel in Fourier space as a power series: 

where c is a completely symmetric tensor. Now we consider S[<p'] with integration variable x', then change to x' — > x 
where x' = x + a with position dependent a. The corresponding Jacobian is dx'^/dx 1 ' = 5^ + d^a v + 0(a 2 ), and the 
transformation of the derivative reads d'=d^ — (d fl a v )d u + 0(a 2 ). We obtain 



(17) 



= E c (*n) /^C 1 + 9a)if(x) ({8^ - d^a^Hd^) . . . (id^ n - d^a Vn id Vn ) <p(x). 



We need the difference linear in da. We find 



SS = d 4 xdaC - E / d 4 x(p(x) E^i) ■ ■ ■ {id^-i^^Hd^id^^) . . . («? M J </?(x). 

^ i=l 

Now we perform partial integration in the second term, and find 

„ n 



Since 



is completely symmetric in the indexes, we can put i — > n and find 



n-1 



J i=l 
Introducing the symmetrized derivative notion we find 

-E™"" 1 " / ^aV(x)[( i 9, 1 )...(za AI „_ 1 )] sym R^(x)) = / ^o"^) 



1 9/C(p) 



Finally 



so we find 



5S = ~ / d i xd"a u {x) { g^ip{x)lC{id)ip{ x) — ip(x) 



dK{p) 



p^-id 



2 d PfJ 



{id v ip{x)) 



p—Hd 



(18) 



(19) 



(20) 



(id„(p(x)). 



J sym 



T^(x) = -<p{x) 



I 


~dK(p) 






dp* 1 


p^id 



(21) 
(22) 

(23) 



This is the equation we had to prove. 
QED 
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B. Finite temperature expectation value 

Now we will compute the expectation value of the energy-momentum tensor in equilibrium. For that we change 
into Fourier space: 

T ^ {k) = W? {27r)4S{k ~ P ~ qMq) D ^,P) <f(p), (24) 

where the symmetrization of the derivatives is inherited to the momentum dependent kernel as 

n 

\p n ]sy m ^— l Y,^) i P n ~ i - ( 25 ) 
i=0 

We can take its expectation value: 

(VW) = \j^f -^^)%k-p-q)D la ,K{q,p) (<p(q)<p(p)) = \ J^D^Kiq = -p,p) *G<(p)(27r)**(A:). 

... _„ -. _,, (26) 

In the expression of K.(q = —p,p), according to (|25[) . the symmetrization of the derivative can be omitted, so we 
simply denote 

D^Kip) = £V ~ 9iu>K. (27) 

Since (|26j) is proportional to S(k), in real space it is position-independent: (T^ u (x)) = (?),„). We remark that if the 
kernel is relativistic invariant (i.e. depends only on p 2 ) then (T^) is a symmetric tensor. 

To determine the propagator at finite temperature we use KMS relation |19j. We note here that the formal 
correspondence between the time translation e~ lHt and statistical operator e~@ H is still valid, and so the KMS 
relation remains true. With the KMS relation we find: 

(T/tv) = 5/7S4 D ^(P) »(Po)<?(P). (28) 

where n(po) — (eP Va — 1) _1 is the Bose-Einstein distribution. Using the fact that K is symmetric real function, D^K, 
is also symmetric real function, while g(po) is a real antisymmetric function. Therefore we can average the positive 
and negative energy parts as 

n(po)g(Po,P) ( ^+n(po) ) g{po,p), (29) 



2 

which is a completely symmetric function. We can then restrict ourselves to the positive frequency part as 

+ . + 00 „ 

dpo f rf 3 P (3Q) 



(T^) = J D^Kip) Q +n(po)) 9(P), where J = J^J 



p 



(2tt)S 



The leading 1/2 represents the vacuum energy. To obtain a finite expression we have to renormalize the expression 
which means "normal ordering" , ie. subtract the expectation value at zero temperature: 



(T^) ren = J D^K.(p) 
p 



+ 

n(po)g{p) + ^Sg{p) 



where 8g(p) = g(p) - g Q (p). (31) 



In a rotationally invariant system only the diagonal elements survive. The energy density is the 00 component of the 
energy-momentum tensor: 



DK{p) 



n(po)g(p) + ^Sg(p) 



DK(p)=p ^-K, (32) 
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This is a completely finite expression for the energy density. This formula is the main result of this paper. 
In case of relativistic invariance the kernel is a function of p 2 alone, then we can write 



DK{p) =2p 2 ^- AC. 



dK 
'dp 2 



(33) 



Note, that (|32[) is non-linear functional of the input spectral function, since apart from the explicit g, the kernel 
K. also depends on it. This nonlinearity has the interesting consequence that the result is independent of the overall 
normalization of the spectral function. Indeed, if we change the normalization of g, ie. we multiply the original 
spectral function by Z, then Gr — > ZGr and /C — > K./Z, and so Z drops out from the result. This fact can be 
interpreted that thermodynamics is not sensitive to the normalization of the fields, only on the density of states. 

Once we know the s(T) relation, we can use thermodynamics to calculate other quantities. For the explicit formulae 
we will not consider thermal variation of the spectrum, and so assume 5g = 0. In that case the pressure is 



+ 

p = -T - DJC(p) In (1 - e-a*) <?(p), 
J Po 



(34) 



and the entropy density reads 



T 

s = [ — DK{p) pp n(p ) - In (1 - e~^) 
J Po 



q(p)- 



(35) 



These expressions are also well defined finite quantities. 



C. Number of degrees of freedom 



The real observables are the energy density and the corresponding thermodynamical quantities. In order to help 
understanding the physics of the effective system, other, more loosely defined quantities can be useful like "number 
ofdof". 

For stable free particles this notion is well defined: it is the number of energy levels at fixed quantum numbers, 
in particular at fixed momentum. The spectral function q(oj, k) consists of Dirac-delta peaks in this case, and the 
number of dof is the number of the Dirac-delta peaks. The heights do not matter, they come from the normalization 
of the operator with which we measure the spectrum. Note however, that the normalization independence means that 
the number of dof must be a nonlinear functional of the spectral function! 

A general input spectral function can have arbitrary structure, usually broadened peaks and other "continuum" 
contributions. The question of the number of dof is thus not well-defined there. We may hope, however, that plausible 
definitions yield approximately the same physical picture. 

A plausible choice is, for example, to compute the energy density e of the system. For a free theory with A^of 
degrees of freedom the energy density is Ndof£o(jn, T), where £o(to, T) is the free energy density of a gas with mass m 
at temperature T . This quantity, therefore, is proportional to the number of dof, so we can define a "thermodynamical 
dof" at any temperature. The definition depends on the mass m, unless we go to asymptotic high temperatures, where 
N doi (T -> oo) = 30e(T -> oo)/(7r 2 T 4 ). 

Another proposal is to use the formula of the energy density and use the ideas of Williams and Weizsacker equivalent 
photon number [22j]. The corresponding definition is the "zero temperature equivalent dof" with the Ansatz: 

oo 

N dol = [^-DK(p)g(p), (36) 
o 

This is a dimensionless, temperature independent quantity which, as we will see later, in case of the discrete spectrum 
yields indeed the number of energy levels. 



V. INDEPENDENCE OF THE DEGREES OF FREEDOM 



We have all the formulae to answer physical questions. The main issue will be to compare the results coming from 
the independent resonance approximation and the one-field representation. 
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A. Dirac-delta peaks 



The first problem that we should check is whether a system that contains N stable excitations with different 
dispersion relations can be represented with N free fields. According to fl7| we expect that it must be true. 

Consider now a spectral function consisting of ./V Dirac-delta peaks with dispersion relations Wj(p) (wj ^ loj for 
i 7^ j), allowing also different normalization: 

g(p) =J2^i5(p - Ui(p)). (37) 



In the independent resonance approximation we assign the following Lagrangian to this spectral function: 

2' ; 



N 1 

Cindep = 2J yPiiido - U)i(id))(pi. (38) 



The normalization is irrelevant, since we can rescale the fields to achieve the standard form. Using this Lagrangian 
and measuring the spectral function of the operator ip = J^. \fZltpi we indeed recover (|37| . 

This Lagrangian represents an iV-dof system. If we calculate the thermodynamical observables, the energy density 
and the pressure we obtain sum of the partial energy densities and pressure: 



N N 

e 

«=i 



where 



P ^i(2^3^ n(Wl) - (40) 

Let us check whether we obtain the same result in the single-field representation, where - correctly - we do not 
introduce new conserved charges into the system. Now we have a single field; from the Kramers-Kronig relation (|10p 
we find for the retarded Green's function 

G * = E TXT- ( 41 ) 

^ Po - uj.i + is 

The kernel is its inverse: it is now a complicated, nonlocal object which have zeroes at po = Wj(p) dispersion relations. 
Therefore we immediately see that JCq = 0. The nonzero contribution to the energy density (1321) . pressure (f34f or 
number of dof (I36|) is the derivative of the kernel multiplied by the spectral function. In case of the energy density 
we need 

dIC po dG R pq ^ Zj 

PQ dpo ~ g\ dpo ~ g\ 4- (po - c^) 2 ' [ ' 



The spectral function evaluates this formula at each po = Wj. We find 
hm po-^— = hm 



dpo po^w, Zf/(po — uJi) 2 + finite 
That means 



(po — UJi) 2 + finite 



Zi 



(43) 



N / Qtr\ N 

DK(p) Q(p) = V lim [p — ) Z i (2TT)S{po-uj i ) = Y / uj i (2n)6(p -uj i ). (44) 

The wave function renormalization factors disappear, and the energy becomes u>i. Therefore the number of degrees 
of freedom from (f36|) is reduced to the integral 

N 00 

iVdot = V / ^ - Wi(27r)*(po - Wi) = JV, (45) 
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supporting the N independent degrees of freedom approximation. The formula of the energy density (|32|) yields 

N } N f d 3 p 

£ = ^2 j n (Po)wi(27r)5(p - cJi) =22 J J2^j3 u}in ( u; i), (46) 



1=1 p 



which agrees with (|39[) with (|40[) . Similarly for the pressure, using Q34p. we recover the pressure formula in (|39p with 

una. 

We can conclude that for a spectral function which has N Delta-peaks, the independent resonance approximation 
exactly reproduces the results of the one-field representation. This also agrees with the earlier results 



B. Single particle with finite width 



Now let us start to study the problem of overlapping peaks. The basic building block is the case of a single particle 
with finite width. There the spectral function has the (relativistic) Breit-Wigner form 

S(JP) = T2 r2 4P °2\2 i A 2r2 " ( 47 ) 

From this spectral function we obtain the retarded Greens function and the kernel as 

Gr(p) = , - * j ' =pg-r 2 -w», D1C=pI+T 2 +lj 2 p . (48) 

The formula for the number of degrees of freedom ([36]) can be exactly evaluated: 

iVdof -i 2^ po (^-r 2 -^) 2 + 4 p 2 r 2 i; ^ 



it is independent on the width of the Lorentzian. 
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FIG. 1. The T 4 normalized (a.) energy density and (b.) pressure of a system with spectral function of a single Lorentzian, 
compared to the same mass, zero width case. SB means Stefan-Boltzmann limit. 



For the energy density and pressure we have to do numerical calculations. The result can be seen on Fig. [TJ For 
comparison we have depicted the zero width (ie. Dirac-delta) case. We can see that despite the relatively large with 
(r/m = 0.2), the energy density and pressure are very close to the free case, especially in the low energy region. This 
makes understandable that the hadron resonance gas approximation works well in the low energy regime even for 
particles with relatively large width like the rho-meson. 
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C. Two Lorentzians 

Now let us take two Lorentzians in the same quantum channel. The spectral lunction which we use consists of two 
independent Breit-Wigner functions: 



q(p) = 



4:Z 1 p T 1 



4z 2 p r 2 



(p§-r?- w ?)* + 4p8r? 



(50) 



The spectral function is shown on Fig. [5] Now the spectral function can be exactly represented by two finite width 




0.5 1 1.5 



2.5 3 



FIG. 2. Spectral function consisting two Breit-Wigner form with the same width (F = 0.2) and with masses m = 1 and 2. The 
dashed lines show the two Breit-Wigner ingredients of the spectral function. 



particles. The independent resonance approximation then would yield the sum of the single Lorentzian results: the 
number of degrees of freedom is 2, and the energy- density and pressure is the sum of the one-Lorentzian results: 

E = £l-Lorentzian('mi) + £l-Lorentzian(m2) and p = Pi- Lorentzian (™1 ) + Pl-Lorentzian('m>2) ■ 

Now we work out the result coming from the one-field representation which does not introduce new conserved 
quantities. The corresponding retarded Greens function is: 



G R (p) 



(pa + Wi) 2 - uj 2 (po + iT 2 ) 2 



,2 : 



(51) 



its inverse is the kernel of the Lagrangian, according to (1111) . From there we can calculate the number of dof (I3GI) . 
the energy density (1521 and pressure (fM]) . These letter need numerical calculations, the number of dof, however, can 
be worked out analytically, and it reflects well, what is going on when the two peaks start to merge. 

By representing the spectral function with its pole structure, the formula for the number of dof p6[) reads in this 

case iV dof = A^ ] f + N$$ where 



N. 



(i) 

dof 



dp DK(p) Z x 
2ir po Aimi 



1 



1 



Po — iTi — mi po - iTi + mi po + iVi — mi p + iTi + mi 



(52) 



(2) 

and NiJt is a similar expression with 1^2 index change. We can close the contour of the integration from above (or 
below, the result is the same) to find 



N {1} - 

Jv dof — 



Ami 



DK[p) 



Po 



DK(p) 



p =m 1 +iT 1 



Pa 



-mi +1^1 



(53) 



Since DK.(—po) = DK,{j>o), the above expression is real. The actual expression is rather complicated, but with help 
of Mathematica one can work out the analytical formulae. The generic formula is rather lengthy, here we analyze the 
Z\ = Z2 and Ti = T 2 case. Then the result reads 



N, 



(i) 

dof 



(2048r 8 + {mf - m\f + 16r 2 (m? - m\) 2 (m\ + ml) + 512r 6 (5m 2 + 3mf) + 32r 4 (23mf + Mm\m\ + 7m%) 



(64r 4 + (ml - ml) 2 + 16r 2 (3m 2 + m 2 )) 2 



(54) 
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FIG. 3. The number of degrees of freedom at fixed masses but changing width of the two Lorentzian peaks. The solid line is 
the definition, the dashed line comes from the approximate formula (|55ll . 



(2) 

and for N^ of we obtain the same formula with 1h2 change. When we choose mi — 1, m 2 — 2 then we can plot the 
degrees of freedom, it is shown in Fig. |3] To understand better these formulae it is worth to work out the approximate 
formula valid in the T, \m\ — m 2 \ <C TOi,2- We obtain 



1 - y 2 4r 

N do{ = l + — — where y=- f . (55) 

(l + jrr \m 2 -mi\ 



Although it is just an approximate formula, it approximates well the correct expression, as Fig. [3] shows. 

The robust feature of the number of dof is that for small width it starts at 2: this is the case of two Dirac deltas 
analyzed in Subsection IV Al On the other hand as Y — > oo the number of dof goes to one. The two peaks merge, 
and there is no way to say whether they were two separate peaks or not. Then the situation is similar to the one 
we considered in Section [II] We can therefore trace how a two-degrees-of-freedom system becomes a one-degree of 
freedom system dynamically, without changing the number of representing fields manually. 

At non-extremal cases, ie. when < T < oo, the number of dof is not a well-defined physical observable. As we 
can see, with the present definition (1361) the number of dof can be smaller than one. Still, in order to determine the 
width when the the originally two-dof system turns into a one-dof system, this definition gives probably a sensible 
answer. Therefore we can claim that already at Y / Am « 0.2 we have an approximately one-dof system, although in 
the spectral function the two peaks are clearly identifiable yet (cf. Fig. [2]). 

We remark here that this phenomenon is closely related to the resolution of the Gibbs paradox. The paradox is 
that if in a two dof system we continuously make the "difference" of the dof to disappear, then for arbitrarily small 
difference we observe a 2-dof system which suddenly, non-analytically changes to a 1-dof system at zero difference. 
The resolution, according to the above analysis is that this change is in fact analytical if we consider finite lifetime 
particles. A similar conclusion is drawn in 23]. Physically it means that we cannot perform a sufficiently long 
measurement which is accurate enough to identify the two separate energy levels in case of finite lifetime. 

Although the number of dof describes well, what is going on, but we still have to see, how it is realized in the 
physically observable quantities. To this end we determined the energy density and pressure for the case of small 
and large width. The result can be seen on Fig. 01 The exact result for Y = 0.2 is the curve labeled by Hi.). If 
we assumed that the two Lorentzians are independently contribute to the thermodynamics, then we would obtain 
curve ii.). For comparison we also show the results of the small width case (ie. two Dirac-deltas) as curve i.) and 
the result of one Lorentzian case with mass m — 1.2 and width Y = 0.2 as curve iv.). It can be clearly seen that in 
the overlapping Lorentzian case the independent resonance approximation fails. While the results of the independent 
resonance approximation is close to the two Dirac-delta peaks case, the exact result is closer to the one Lorentzian case 
with appropriate parameters: ie. a one-dof approximation. Therefore the results of the thermodynamical calculation 
which provides physical observables supports the conclusion coming from the analysis of the number of dof (|36[) : for 
r/Am = 0.2 the system with two Lorentzian peaks practically behave as a one-dof system. 



D. Thresholds 



We have seen in the previous subsection that two Breit-Wigner resonances in general are combined in a non-linear 
way. This phenomenon is in accordance with the scattering theory predictions |15| , where the S- matrix contributions 
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0.5 1 1.5 2 2.5 3 0.5 1 1.5 2 2.5 3 

T/mj T/mj 

a.) b.) 

FIG. 4. The T 4 normalized ( a.) energy density and (b.) pressure of a system with spectral function consisting of two Lorentzian 
peaks with mi = 1, vni — 2 and common T width, cf. eq. (|50(l . The result for T = 0.2 is in both figures denoted by Hi,), the 
other curves are here for comparison, i.) is the result of the two Dirac-delta case, ii.) is the result of the independent resonance 
approximation, iv.) is the result of one Lorentzian case with mass m = 1.2 and width V = 0.2. SB means Stefan-Boltzmann 
limit, 2SB is twice the Stefan-Boltzmann limit. 



of the resonances are non-independent, but they satisfy in general a complicated unitarity constraint equation. We 
should go further and study, what happens when the spectral function can be represented by several resonances. Does 
this nonlinear reduction of the number of dof take place also in this case? 

The generic formulae for the energy density, pressure and number of dof are shown in (j32]) , (|34j) and (|36|) , they can 
be evaluated for an arbitrary spectral function. Now let us choose a simple threshold spectral function 



Q(p) = 8(p 2 -M 2 )^--^. (56) 
The spectral function is shown in Fig. [5] In principle one could represent this spectral function, with some prescribed 




FIG. 5. The spectral function of a single threshold 



accuracy, with a number of massively overlapping Lorentzian peaks. Since the overall normalization drops out from 
the formulae, we do not even have to care about the correct normalization. The threshold is, however, very far from a 
particle-like spectral function: there is no peak, it does not vanish at infinity. We actually do not even expect that it 
describes a physical system on its own: rather it serves as a mathematical example for the strongly correlated peaks. 

When we evaluate the number of dof using (f36|) we obtain for m = 1 the value Ndof = 0.98. However, when we use 
the formula for the energy density and pressure we obtain the curves of Fig. |6l These curves are very far from the 
multi-dof expectations. The most striking difference is in the amplitude: instead of having several times the Stefan- 
Boltzmann limit, we are below the 0.01 x SB limit. This further supports the hadron resonance gas approximation 
where one neglects the contribution of the branch cuts to the thermodynamics. We see that in this pathological case 
the number of dof coming from (|36|) was also not informative, this seems to work only for particle-like systems. 

If we want to describe what is going on using the language of multiple resonances, we can say that the contributions 
of the single resonances interfere with each other so strongly, that their net contribution to thermodynamical quantities 
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FIG. 6. The energy density and pressure of the pure threshold spectral function. SB/100 means the hundredth of the Stefan- 
Boltzmann limit. 



is almost completely washed out. To understand how is it possible we recall that the single resonance contributions 
are complex (cf. |15|). but their sum finally must give a real number. So we have a large number of uncorrelated 
complex phases, and study their sum: it is then plausible that finally we obtain a near-zero net contribution. 

We emphasize that this last example is non-physical. It serves for providing a warning signal that it can be very 
misleading if one tries to represent a spectral function with a large number of heavily overlapping peaks. 



VI. CONCLUSIONS 



In this paper we discussed the question of what kind of effective model can be built on an experimentally determined 
spectral function. The spectral function in general can consist of multiple finite width peaks and non-particle-like 
continuum contributions. One procedure can be to represent the spectral function with the sum of several Brcit- 
Wigner spectra, and assign a new dof for each of these Breit-Wigner peaks which interact with each other in a way 
that the correct widths come out (independent resonance approximation). However, in this case each of the new 
degrees of freedom at tree level introduces a new conserved quantity (the corresponding particle number), which may 
alter the dynamical behavior of the system. While for large separation this procedure is known to work (cf. [17|). 
with massively overlapping peaks the independence of these resonances is questionable. 

A more conservative way is to find a representation which does not use extra, non-physical dof. As a price we 
have to describe multiple peaks with a single field. This can be done consistently, defining carefully a non-local 
effective theory, and checking all the consistency requirements. In this case the number of dof is not the number of 
the representing fields (which is one), we can define it from the spectral function. 

Having a physical representation at hand, we can study how well the independent resonance approximation works. 
We have found that for infinite lifetime excitations (ie. a spectral function which consists of discrete Dirac-delta 
peaks), the two methods provide identical results. So, as we have expected, infinitely well separated peaks indeed 
describe independent dof. But if the peaks are finite width Breit-Wigner resonances and the peak separation becomes 
comparable with the width, this nice agreement starts to fail. As we have shown, the relevant parameter is T/ Am: a 
two-peak spectral function with T/ Am = 0.2 already behaves from the point of view of thermodynamics as a one-dof 
system. An even more striking example for the non-independence of the overlapping Breit-Wigner peaks is the case 
of a mathematical construction of a threshold-like spectral function. Although it could be represented with a sum of 
multiple, overlapping Breit-Wigner resonances, the energy density of this model is still about 100-times smaller than 
that of a one-particle system. 

We may conclude finally that a multi-dof representation of the spectral functions can yield very misleading results. 
According to the numerical evidences of this paper, the peaks of the spectral function represent independent dof only 
in the case when the width of the peaks is smaller than about 0.1 x the separation of the peaks. In other cases one has 
to use a one-dof representation in order not to change the original symmetries of the system, and obtain physically 
correct results. 
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